Generalized Elastic Model yields Fractional Langevin Equation 
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Starting from a generalized elastic model which accounts for the stochastic motion of several 
physical systems such as membranes, (semi)flexible polymers and fluctuating interfaces among oth- 
ers, we derive the fractional Langevin equation (FLE) for a probe particle in such systems, in the 
case of thermal initial conditions. We show that this FLE is the only one fulfilling the fluctuation- 
dissipation (FD) relation within a new family of fractional Brownian motion (FBM) equations. The 
FLE for the time-dependent fluctuations of the donor-acceptor distance in a protein, is shown to 
be recovered. When the system starts from non-thermal conditions, the corresponding FLE, which 
does not fulfill FD relation, is derived. 
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PACS numbers: 05.40.-a, 02.50.Ey 

Introduction.- Continuum elastic models have been 
extensively used in statistical mechanics to study the 
dynamics of real physical systems. Examples include: 
(semi)flexible polymers [ll-Q, membranes p, JJ-Q, 
growing interfaces |7 Hlll, fl uctuating surfaces [l^ and 
diffusion- noise systems [13|. In this paper we consider 
the following Markovian equation for a generalization of 
the accounted elastic models {generalized elastic model) 



rh(f',i)+J7(f,t) 



(1) 
for the dynamics of the £)-dimensional stochastic process 
h in the d— dimensional infinite space: the internal d co- 
ordinates are represented by x and the Gaussian white 
noise satisfies the fluctuation-dissipation (FD) relation 

(77, (a-, t) ?7fc (f, t')) = 2kBTK (f - f ) 5j k6{t ~ t'). (2) 



with j, fc e [1,1?; 

fined via its Fourier transform as J'qi -nvm 
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r := -(-A) 



The Riesz fractional operator is de- 
or, in terms of the Laplacian A, as 

/2 



15 



In the following analysis we 
study two classes of hydrodynamic interactions: long 
ranged, A (r) ^ \r\ as f — ;■ cxd, with Q < a < d, and 
local, A (r) — 6"^ {f). In the non-local case, if a = d we 



take A (f) 



-kT 



where a is a microscopic cutoff. 



Long range hydrodynamic interactions. The hydrody- 
namic interactions are often represented by the equilib- 
rium average of the Oseen tensor, which in an embedding 
de-dimensional space (de > 3) reads [ij: A (r) ~ |r| ~ ^. 
Examples are: (I) fluid membranes [2, |j-l6|, whose 
height h (x, t) of a point x on the 2-dimensional (planar) 
base surface is moving in time according to ([T]), with 
a — D + d~2 — 1 and z = 4 as derived from the Hel- 
frich bending free energy for small deformations [16|. 
(II) Semiflexible and flexible polymers ' models, where h 
represents the 3 spatial coordinates of a polymeric seg- 
ment (bead), while x is the strand's 1-dimensional inter- 
nal coordinate (curvilinear abscissa). For semiflexible fil- 
aments [2I, y the bending elastic energy associated with 
the chain's deformation implies 2 = 4 17[ and from the 



Oseen tensor formula we get a = 13 — 2 = 1; for the 
fiexible polymers, often referred to as Zimm model [l|, 
the free energy contribution solely comes from the elastic 
term, i.e. z = 2, and a = 1/2 in solvent. 

Local hydrodynamic interactions. Examples of the 
case where hydrodynamic interactions are completely 
screened out are: (I) the Rouse model equation for poly- 
mer dynamics [1|, once one sets D — i, d = \ and z = 2. 
(II) Single file system: recently it has been shown [l8| 
that the dynamics of a gas of Brownian hard rods on 
a line can be mapped onto the harmonic chain problem 
(z = 2), where h(x, t) stands for the position of the x-th 
particle on the 1-dimensional substrate at time t. (Ill) 
Fluctuating interfaces [7, ^] , where h plays the role of 
a scalar field (mostly the height of a rough surface in 
d dimension) which is subjected to a non-standard elas- 
tic force embodied by the fractional derivative of order 
z. This is actually the generalization of the Edwards- 
Wilkinson equation for the fluctuating profile of a granu- 
lar surface, for which d = 2, z = 2 [sf. In systems such as 
crack propagation [lO[ and contact line of a liquid menis- 
cus [llj d = 1 and the restoring forces are characterized 
by z = 1. If instead h is meant to be a step, namely a 
line boundary at which the surface changes height by one 
or more atomic units, the value of z in eq.([T|) is found to 
be z = 2, 3 or 4 (d = 1) according to the character of 
the atomic diffusion [l2[. (IV) Diffusion-noise equation 
[13J : in this case h represents the density field on a d- 
dimensional surface x and z = 2. 

The aim of this manuscript is to derive the non- 
Markovian time-FLE for the field hj at a given position 
X (probe particle), given that the whole system's dynam- 
ics obeys the space fractional Markovian equation ([1]). 
We interpret the FLE as a particular case of a broader 
class of stochastic equations for FBM, i.e. the generalized 
fractional Langevin equation. 

Autocorrelation function.- We hereby calculate 
the /i-autocorrelation function (6h(x,t) dh(x,t')) = 



[hj (x, t) - hj (x, 0)] [hj (x, t') - hj (x, 0)]). 

in space and 

the solution of 



ing the Fourier 

hj (q, uj) — J d'^x dt e 



transform 

*(«-^""*'/ij(f,t). 



Defin- 
time as 




FIG. 1: (Color online) Schematic picture of the GFLE (/i, H) 
parameter space. 



©is 



hj {q, w) = 



V] (g. ^ 



-ILJ + A{q)\q]' 



(3) 



where A (q) stands for the d-dimensional Fourier trans- 
form of the hydrodynamic interaction term, reading 

Al-Z) - '%r '%f,i'' \r-' = A\qr~'; a a ^ din 
first approximation we can neglect the logarithmic cor- 
rections in the hydrodynamic term's Fourier transform 
[31, i-c. A (g) ^ const. In the local hydrodynamic situa- 
tion A{q) = 1, which corresponds to put A = 1 and a — d 
in the long-ranged hydrodynamic expression: this substi- 
tution (which is not to be intended as a limit) allows us 
to easily shift from power-law to local hydrodynamics 
throughout the following analysis. Moreover, due to the 
isotropy of the problem under study, we can drop the 
label j in ([3]). The ft,- autocorrelation function is readily 
obtained: the general expression looks like 



{Sh{x,t)Sh{x,t')) ^K 



t^+t'^ 



|t-tf 



(4) 



where the anomlaous diffusion exponent /3 and the diffu- 
sion constant K are defined as 




a+z — d 
(27r)''r(d/2) z-d 



(5) 



Henceforth we will concentrate on the case z > d, which 
does not need any regularization of the (f integrals in the 
inverse Fourier transform. The cases z < d and z = d will 
be reported elsewhere. From (j4]) it is apparent that the 
probe particle placed at a given x performs a fractional 
Brownian motion (FBM) which is always subdifFusive, 
namely {5'^h{x,t)) ^2KtP. 

Fractional Langevin Equation.- We now proceed to 
derive the fractional Langevin equation (FLE) for a 



probe particle placed at position x, which fulfills the 
FD relation. In what follows we show that this is the 
only fractional stochastic equation physically relevant 
for the considered system ([T|). We first include long- 
ranged hydrodynamic interactions. We multiply both 
sides of the Fourier solution ^ by A'+(— iw)^, where 
We then define the following function 






~P)- 



$ {x,uj) 



d'^q e-'9-^ {-iLu)^ 



Inverting the Fourier transforms in space gives 

{^iujfK+h{x,uj) =C(^,w) 



(6) 



(7) 



where we have introduced the fractional Gaussian noise 
(fGN) Cix,uj) = K+ Jd'^x'r]{x-x',uj)<!>{x',uj). In 
time domain the previous equation takes the final form 
of a FLE [i,[ll[li[2i,[23]: 



K+Dlh{x,t) =C{x,t) 



(8) 



where D^ is the Caputo derivative defined as [15|, |2C 



Difit) = 



1 



r(i-/3) 



dt^j^ifitro<p<i. 



(9) 
The noise in ([5]) can be shown to satisfy the FD relation. 



{({x,t)C{S,t')) = kBT 



X^ 



T{l-(3)\t-t'\ 



/3- 



(10) 



In the absence of any hydrodynamic interaction the 
anomalous diffusion exponent /3 in (jH]) and ([T0|) is /? = 
1 — d/z [8|. Along these lines, the FLE for the case 
d — 1, z = 2, A{x — x') = 5{x — x') has been obtained 
in [l-S]: our derivation can thus be viewed as a general- 
ization of the technique used in such model. We empha- 
size that the spatial correlations appearing in the model 
([1]) are translated into time correlations described by the 
fractional derivative ([9|), together with the space-time 
correlations of the noise C [x, t). As a consequence, any h- 
correlation function can be calculated starting from ([8|), 

e.g. ( h{x,t)~h{x',Q\ ) has been studied for fluc- 
tuating membranes [6|. We note that the underlying 
assumption in the expression ([8]) is that the system has 
reached the thermal equilibrium at i = and the sys- 
tem's configuration is drawn from the stationary Gibb- 

sian probability distribution '^ e ^ V 3i»=i ' 

na. 

Now, it is possible to recast the result (fTOj) in the same 
fashion as the fGn correlation function of a FBM 



22 



tf ^-^ with H 



1 



^. It stems 



from ((H) and ([5]) that the /i-autocorrelation function can 
be expressed as {5h{x,t) Sh{x,t')) ex i^-an ^ ^/2-2H _ 
\t — t'\ ~ , which is at odd with the corresponding stan- 
dard FBM quantity [22|, whose the exponent is 2H. 
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FIG. 2: (Color online) First passage time (FPT) distribution 
of the tagged single file particle initially placed at distance L 
from the adsorbing boundary. The system consists of 2 x 10* 
Brownian particles subjected to hard-core interaction. Simu- 
lation' details can be found in 18. , 26] . The parameters are: 
the file's particle density p — 0.25, the temperature fcsT — 1 
and the damping 7. The theoretical prediction r^ t~^'^^ has 
been drawn for the reader convenience (dashed line). 



However this is not surprising, since the two processes 
obey two different stochastic fractional equations: i) the 
FLE (|8]) for systems which fulfiU the FD relation, and ii) 
the usual equation 



dBHJt) 
dt 



Cnit) 



(11) 



forFBM ^2^. 

Generalized Fractional Langevin Equation.- Let us 
now generalize eq.(j8]) and (fTTj). Consider a stochastic 
process G{t) governed by the following dynamical equa- 
tion 



DlG{t) - Ch (t) 



(12) 



where (h {t) is a fGn which satisfies {C,h (t)) — and 
{Ch (t) Ch {t')) ~ C |t - t'f"~^ for |i -t'\^oo with < 
H <l, H ^ 1/2; C<OifO<i7< 1/2, and C > ff 
1/2 < H < 1. For H = 1/2 (C > 0) Ch {t) is the white 
Gaussian noise [23|. The fractional derivative has been 
defined in ([5]): it is immediate to recover equation ((TT]) 
in the limiting case /i = 1, once that D\_ = d/dt and 
G = 2H{2H - 1), and also the FLE §\), setting K+ = 1, 
for systems satisfying the FD relation: n — 2 ~ 2H and 
G = kBT/T{2H — 1). The autocorrelation function can 
be calculated to yield 



{SG {t) SG {t')) = 



Gsm{TrH)T{2H -1) 



^ ^2iH+t,-l) 



+ t 



sin (vr {H - 

/2(H+fj.-l} 



l))r{2H + 2^i-l) 



-\t-t 



,|2(H+^-l) 



(13) 



Expression P^ shows that G{t) is a FBM with Hurst 
exponent Hfbm = i? + /i-l {1 < H + fi< 2 ^). As 



a consequence, systems for which H + n < 3/2 exhibit 
subdiffusive motion, which instead is superdiffusive for 
H + fi > 3/2. It is interesting to note that the class of 
FBM for which the FD relation holds can be only subd- 
iffusive. These results are summarized in Fig. [Tj . 
An important corollary of ([T^ states that any statisti- 
cal property which is shown to be valid for a given pair 
of values /i' and H' , is automatically valid for any other 
pair (^", H") which satisfies ^i' + H' = /t" + H". The 
demonstration is straightforward: it is sufficient to note 
that, since G{t) is a Gaussian process, it is fully spec- 
ified by the correlation function (|13l) . As an example, 
take the first passage time distribution (FPT) in a semi- 
infinite domain for a FBM, which decays asymptotically 
like r^ t""'""-^ (1, [ii]. We immediately get that the 
FPT distribution for a process which is solution of p^ 
is given by ~ i'^+^~'^, which in turn reads ^ t^~'^ for sys- 
tems obeying to ([1]). We numerically support this result 
as shown in Fig. [5]. The FPT distribution of a tagged 
particle in single file system, which has been shown to 
be described by a FLE with i7 = 3/4 [H, [l^l, since 
/i = 2 — 27J, attains the ~ t~^-^^ asymptotic behavior. 
On the other hand, given an FBM process with 



{5'G{t)) 



2Hfe 



there is no chance to determine 



the correct pair (//, H) among the GFLEs p^ for which 
fi+H—1 — HpBM- Nevertheless, introducing an external 
potential does the trick: for instance, adding a constant 
force F on the right-hand side of ((T2|) gives {G{t)) = 



F 



t" 



r(M+i) 
the relaxation 



while an harmonic 



force —i}'^G{t) leads to 
= (G(t)G(0))/(G2(0)) ^ 



E„.A \~ (t/tn)''!, in case of deterministic and thermal ini 



n-^/f" and E. 



'JAM 



^^ {G{t))/G{0) 

tial conditions respectively, where to 
stands for the Mittag-Leffler function [20|. 
One might question the uniqueness of the FLE ^ among 
the whole family of GFLEs for the probe particle h (x, t) . 
It is possible to show that eq.([8|) is the unique GFLE 
for a process h {x,t) whose dynamics is ruled by eq.([T]). 
The demonstration deals with the introduction of a local 
constant force field Fj5 [x — x*) 6{t) on the right-hand 
side of eq.iH]). Since the system fulfills the FD relation 
([2]), the connection between the average drift of h {x* ,t) 
and its mean square displacement in the absence of force 
is given by the Einstein relation 



{hix*,t))F^F 



{5'^h{x\t)) 
2kBT 



(14) 



The only GFLE which reproduces p^ is that which pre- 
serves the FD relation, i.e. the FLE ([5]). 

Let's now briefiy discuss a practical example of the use- 
fulness of the framework developed here. In Refs 27ll28| 
Sunny Xie and coworkers succeeded in modeling the 
motion of the donor-acceptor (D-A) distance within a 
protein, as the coordinate of a fictitious particle dif- 
fusing in an harmonic potential according to a FLE 
with fractional derivative of order 1/2. In the spirit of 
Refs [29|, l30| , we consider an idealized Rouse chain as a 
model for the protein conformational dynamics. There- 



fore we take D = 3, d— 1, z — 2 and A{x — x') = 
6 {x — x') in (IlJ. The D-A distance vector can be ex- 
pressed as A£)_^(i) = h{xA,t) — h{xo,t), and its 
correlation function by [30| {A.ij_A{t) ■ A£)_^(i')) — 
3{AD-A{t)A£)-A{t')) due to the isotropy of the sys - 
tem. Hence, we can employ the result of Ref. |18| . 
which shows that the generalized Langevin equation for 
the single component A£)_a(0 is 1/2 i?+ A£)_^(t) — 
-wg [Ao-Ait] - Ad-a(O)] + Ci3-A(i) in the long time 
limit , with wq k (xa — xo)^^''^ and Cd-aC^) satisfy- 
ing the FD relation. However we point out here again 
that (A£)_^(i)A£)_yi(i')) can be evaluated directly from 
eq.®. 

Non-thermal initial conditions. -LeVs now assume that 
the initial conditions for the system in ([1]) are given by 



h {x, 0) = 



(15) 



without loss of generality. For systems such as fluctuat- 
ing interfaces |7H9l Il2| or membranes [2,|j46|, eg. p^ 
assumes the interface to be flat at i = 0. In the case 
of a polymer, we can imagine eq. (|15p to be valid only 
for the j-th component, achieving an initial configura- 
tion which is randomly arranged within the plane j = 0. 
For single file systems eq. ([T5|) consists of taking particles 
equally spaced at i = 0. The /i-autocorrelation function 
can be obtained in the same fashion as in the case of 
thermal initial conditions by using Laplace transform in 
time instead of the Fourier transform. A straightforward 
calculation yields 



{6h {x, t) Sh {x, t')) ^K {t + t'f -\t-t'\ 



(16) 



where the value of K and /3 get the same expressions 

as in ([5]). For local hydrodynamics eq. p^ matches the 

result previously obtained by Krug et al. for fluctuating 

interfaces ^] . 

It is easy to show that the FLE expression ^ is still valid, 

with the Ca put o derivative having its lower terminal at 

t 



m, 31 



When attempting to recover the FD 
relation, however, one gets the following form of the noise 
correlation function 



{c{x,t)ax,t')) 



kBTK+ 

r(i-/3) 



|i-<'r -(i + i')" 



(17) 

Expression (ITTI) clearly shows that the noise C attains the 
stationarity ([T0| in the limit (i,i' -^ oo). 

Discussion. - In this Letter we presented the derivation 
of the FLE for a wide class of phenomena, whose stochas- 
tic dynamics is ruled by the generalized elastic model ([I} . 
The introduced framework offers theoretical and practi- 
cal advantages. On one hand, different physical systems 
can be deflned on the basis of a unique index: the frac- 
tional derivative order [universality class). On the other 
hand, the FLE allows to achieve the relevant statistical 
observable by simply solve/simulate a non-Markovian lin- 
ear equation for the probe particle. Finally, from an ex- 
perimental perspective, the FLE description allows the 



straightforward detection of the microscopical parame- 
ters characterizing the system ([ij . 
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